Observability of glueball spectrum in QCD and the width of a resonance 
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We prove a theorem in QCD stating that in the limit of strong coupling, g oo, the observed 
spectrum of glueballs in QCD is the same of a pure Yang-Mills theory, being mixing effects due 
to the next-to-leading order. A full effective theory for QCD is obtained and the width of the a 
resonance decay is straightforwardly computed. This appears as the lowest glueball state. Vacuum 
gluon condensate is computed that consistently support studies on the identification of this meson 
as a glueball. 
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The question of the light scalar meson spectrum is an open problem whose solution would imply a great improve- 
ment in our ability to manage low-energy QCD. This understanding is currently lacking and we have difficulties to 
comprehend the nature of a number of states that appear in laboratory. One of the main open questions is if glueballs 
can be seen and, if they have been seen, what observed states these are. This lack of knowledge has also implied 
serious difficulties into the identification of molecular or tetraquark states that some researchers assume representing 
the light scalar spectrum [ij. This view is challenged by some recent analysis by Mennessier, Narison, Ochs and 
Minkowski on the a reonance 0, Q giving evidence for its glueball nature. More recently, Kaminski, Mennessier and 
Narison gave a stronger evidence from KK decay Q ■ Indeed, one of the most demanding problem is to know what is 
the structure of the lowest resonance seen so far, also presented as f0(600). Initially its very existence was strongly 
debated. Today, there are precise determinations from experimental data of its mass and width 0, Q and is inserted 
into particle listing 13]. Its observed decays are a — )■ tt+tt^ and ct ^ 77 being the former largely dominant. So, a great 
understanding of QCD would be achieved if we would be able to compute the mass and the width of this particle. 
The aim of this paper is to show how to obtain these results starting from QCD. Our main conclusion will be that this 
particle is indeed a glueball and the ground state of pure Yang-Mills theory. This result relies on an essential way to 
see how quarks and glueball state can mix in QCD. There is no serious proof about and is generally believed that, also 
in the low-energy limit, some mixing must happen. But we will prove a theorem stating that the observed glueball 
states should have almost exactly the properties of the states of a pure Yang-Mills theory and mixing is showed to be 
ineffective. 

The line of research started with the introduction of a gluon condensate [1| strongly support our conclusions. This 
means that there must be a deep connection between our results and those obtained by the aforementioned authors. 
We will show that a vacuum gluon condensate is indeed different from zero and has the expected value in agreement 
with recent analysis [^, This proves that all these studies are well founded and their conclusions correct. 

So, we prove the following theorem that will give the concept of ordering needed for successive computations. 

Theorem 1 (Spectral Theorem) At the leading order of a perturbation series of QCD for the coupling g going to 
infinity the observable spectrum of glueballs in QCD is the same of a pure Yang-Mills theory. 

Proof We take the following actions for the Yang-Mills field [ll| 
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with g the coupling constant and a fixing the gauge. Similarly, /"''^ are SU(N) structure constants and c"^ is the ghost 
field. The action for the fermion fields is 



Sra = J d'^x'y^ q{x) 



i-f-{d + ig^A" ) - mq 



q{x) 
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being q the fiavor index and A" the group generators. So, the generating functional for the quantum field theory is 



Z[r7,7;,e,e-,j] = \{[dq][dq\[dc\[dc][dA\, 
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In order to prove the theorem the choice of the gauge is not relevant and so we take one that makes Yang-Mills 
action simpler. We choose a = 1 so that we can write 
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The reason for this choice will be clear in the following. 

We need to make contact with the phenomenology at low energies. In this case we have a large coupling g making 
ordinary weak perturbation theory not applicable. In order to uncover an asymptotic approximation in this regime, 
we consider the formal limit g — > oo that is dual to the weak perturbation limit 5 — > 0. So, for a generic field 
entering into the generating functional, we ask a formal solution series like 



(5) 
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but, apparently, this series cannot be straightforwardly obtained from the generating functional. An obvious way out 
to this impasse has been proposed in [l2, IJ, 14, 15|. We choose to rescale the time variable as t VNgt being Ng^ 
't Hooft coupling. So, we rewrite the actions above as 
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Ng 



d^x 



1. 



and so, a perturbation expansion will produce a gradient expansion being gradient terms of higher order in the limit 
(7 — > 00. Then, for the quark fields we have 



5„ 



70 • idn 



A"' 

2Vn °J ' 2Vn 



d'^x^q{x) 
q 



7. • 



q{x) 



(7) 



Ng 

We can expand the generating functional as 

Z[7^, fj, e, e, j] = / n [dq] [dq] [dc] [dc] [dA] 
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where we have just kept terms of order O i^vNgj and 0(1). Now we do the key observation that the first term of 

the expansion, in the hmit g ^ oo is that of the pure Yang-MiUs field that is O (^y^g^ and effects due to quark 

fields are just 0(1). This implies that, in the given limit, the spectrum of the pure Yang-Mills field must be observed 
with negligible corrections due to the mixing with quark fields. The reason is that quark currents do not drive the 
Yang-Mills field at the leading order that behaves as it would be free. 
Finally, corrections due to gradients just appear at order O ("y^) • 



Some considerations are in order. The limit g oo implies that the leading order term 
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can be evaluated in the semiclassical approximation. So, if we know a proper classical solution of Yang-Mills equations 
we will be able to get both the spectrum and the propagators of the theory at this order. 

In order to have an understanding of this part of the theory one should care about what are the classical solutions 
to be selected to work with in a quantum field theory. Classical Yang-Mills equations admit both integrable and 
chaotic solutions. Here we just put forward a conjecture, to be proved, that a quantum Geld theory can only exist 
with integrable classical solutions. The reasons to believe this rely mostly on the computational opportunities that 
integrable solutions do provide. Classical Yang-Mills equations admit a lot of integrable solutions and this can be 
easily selected with the so called Smilga's choice and the mapping theorem [l6|. These solutions map Yang-Mills 
theory on a quartic massless scalar field theory. For our aims this is equivalent to 

A^^) = v:.m (10) 

with (j){x) a solution of the equation dfcf) + Ng'^(p^ = and e.g. 

-((0,1, 0,0), (0,0, 1,0), (0,0, 0,1)) (11) 

for SU(2) and 

77^ = ((0, 0, 0, 0), (0, 1, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 1, 1, 0), (0, 0, 1, 0), (0, 0, 0, 1), (0, 0, 0, 1)) (12) 

for SU(3) but the choices are a very large number and increase choosing a larger group. Then we can apply a Lorentz 
transformation and we will get a mapping solution between a scalar field and Yang-Mills theory. For our aims it is 
enough to consider the solution of the scalar field equation removing the gradient part. Then, starting with integrable 
solutions we now compute the Green function at the leading order. We have 



SZ 



Z Sj-{x)Sjb{y) 



1 



m 



exp 



[dA][dc][dc]Al{x)Al{y) 



(13) 



iVNg / d^x 



-doAldoA'"' + dac^d^c^ 



\/N \/ N 



1 
4iV 



ahc j;ars a b 



A'^AlA'^f'A'"' 



+ 



1 



Ng 



4 



Excitation 
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Theoretical 


Error 


a 




1.198140235 




0++ 


3.55(7) 


3.594420705 


1% 


0++* 


5.69(10) 


5.990701175 


5% 



TABLE L Comparison for the 0++ glueball spectrum for SU(3). 



We can apply to this the above mapping theorem using the aforementioned integrable solutions giving 



1 



SZ 



1 



^jPj , [dA][dc][dc]7j;jMx)^{y) X 



(14) 



exp|i%/iV.g(iV2 - 1) / d'^x 

o- ' 



Ng 



with really noteworthy simplifications. It is interesting to see how the ghost field becomes decoupled from the Yang- 
Mills field and so it is that of a free particle. These results were already seen in A Green function can be 
considered also in the limit g ^ oo as 14, 



d^G{t) + Ng^Git) = A^S{t) 
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being A an arbitrary constant having the dimension of energy. This is an integration constant of the theory and 
should be experimentally determined but an higher order theory should be able to compute it. The solution is easily 
written down as 



G{t) = 9{t)A 
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and the glueball spectrum is given through the Fourier series of the Jacobi sn function by 



TO„ = [2n + 1) 
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being K{i) 



1.311028777 an elliptic integral and we have identified a = (^^^^^ ^ A^ as the string tension of the 

/2/3(T5f7(3). This has been verified some years ago with 



theory. This implies that, for the same value of A, 175(7(2) 
lattice computations [l7l |. 

This spectrum, as is and together with the theorem above, makes a strong prediction. We can exploit this if we 
identify a set of pure numbers as 



an = ^ = 2n + 1 — — 
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the first ones being given in Tab. |T] compared to lattice computations [3]. The most striking result is the existence of 
a lower ground state that can be identified with the observed a resonance at about 500 MeV 0]. As said above, the 
very nature of this resonance is hotly debated yet [l3| and is not seen in any lattice computation (i^l • From the above 
theorem in QCD together with our solution of Yang-Mills quantum field theory |16j we can conlude immediately that 
this lowest state in the spectrum is a glueball. Further experimental analyses is needed to clarify this point. Anyhow, 
on the basis of the theorem we have just proved, the full glueball spectrum must be seen experimentally. 

We want to pursue such an identification further, a resonance is known to have a quite large width. So, we estimate 
it evaluating the next-to-leading order correction to the solution presented here. Using again the mapping theorem 
[1^ and the integrable classical solutions of Yang-Mills theory, we can write down a full generating functional for an 
effective theory of infrared QCD. We will have 



Z[v,V,Jc^] = / m Ylmm expUVNg{N' -I) j d^x 
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haviii K in troduced {N'^ — — j^j^ri"'^ ■ Then, we consider the following approximation that holds just in the infrared 
limit 21, 23], this is indeed the leading order of a small time expansion. 



0(x) 
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where the multiplicative factor is there to account for rescaling in time and 

A{x -y)^ 5\x ~ y) [G{t2 - h) + G{h - t2)] 



(20) 
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is the Feynman propagator (see ea. (|16p ). Such a small time approximation turns the functional for the scalar field 
into a Gaussian form 1J|. In order to show this, we compute the generating functional taking 



Ng 
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Ng 



d^yA{x - y)j4,{y) 
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where use is made of the aforementioned small time expansion. Working in this approximation, we can neglect terms 
of order higher than second as these imply higher powers of time variable. So, removing the rescaling in time, we are 
left with the following Gaussian approximation 



Z[v,V,.U] ~ exp <^ -{N^ - 1) / d'^xd''yj^{x)A{x - y)j^{y) \ x 



(24) 
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The remaining integral is Gaussian and we will be able to compute it when we know the quark propagator obtained 
by solving the equation 



ijodo + g—liVt I d'^yA{x - y)j4,{y) 



S[j^,x] = 5^{x) 



to give us finally the following generating functional 
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expji J d'^xd'^y'^fjq{x)S[j^,x - y]riq{y)[ 
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and so our final aim will be to compute this functional of j^. Indeed, this can be straightforwardly obtained as 



A" 



S[j<i,,x-y]^e{t^-ty)S^{x-y)exp<ig—joliVi / dt' / d'^xiA{t' - ty - ta;^, x - y - xi)j^{xi) 



At this stage we do the Nambu- Jona-Lasinio approximation on the gluon propagator 23|, [24 1 

3.76, 



A(a; - y) R 
being a the string tension. So, one has 

S[j,p, x-y] = 0(tx - ty)d^{x - y) exp 



-8\x - y) 



A° 
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being the coupling — G^vji /\/47raI and Gnjl — 3.76^ the Nambu-Jona-Lasinio couphng. 
From this we can evaluate the vertex at a tree level as 
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corresponding to the process a 
Heaviside function 



In order to evaluate this vertex we use the following expression for the 
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We have to go from this n-point function to the probability amplitude introducing pion fields. We do this with the 
following rule: we remove the contributions from the Heaviside functions with the product rriaf^, assuming a at rest 
and rUa its mass and /jr the pion decay constant being about 93 MeV. This corresponds to LSZ reduction. So, when 
we take the trace of the square of the amplitude, we will be left with 



\M.ff^N(N'-l)DGlmlf^ 
that we specialize to iV = 3 and D — A. This will give the rate 



6 G 
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(33) 



that is in agreement with recent derivations from experiments 5] being = 255 ± 10 MeV for a mass ma = 

484 ± 17 MeV when ag ~ 1.4. A similar conclusion can be drawn also with respect to the derivation given in 
presenting a mass ma = 44llg® MeV and Ta/2 — 2791^2 5 MeV, for as ~ 1.8 showing that both results arc consistent 
each other with respect to physical values of the strong coupling constant. We have kept the string tension fixed in 
both the computations at (440 MeV)"^ but this choice may be too tightening with respect to the mass of the particle. 

Our next step is to show how a gluon condensate emerges from our strong coupling computations. This step is an 
essential one as several authors used this concept to show that a meson is indeed a glueball with a mass in agreement 
with the one we obtained above. Indeed, we can show that there is a gluon condensate and its value is in close 
agrement with recent estimations and, in any case, very close to the value estimated by Shifman, Vainshtein 

and Zakharov ^8|. So, we take SU(3) for the gauge group and apply the mapping theorem for Yang-Mills theory and 
a quartic scalar field theory 16] giving 

(G-G) = 16{{d(f>f) ~ 96TTas{ 
This equation implies that we have to compute the following integrals 
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that are not defined unless we introduce a cut-ofF. This is normally done in this kind of computations and, as seen 
above, QCD has a natural cut-off by its own that makes these integrals meaningful. So, for the gluon propagator, we 
introduce also a dependence on the space momentum after resumming the gradient part of the action that can always 
be done and introducing a mass M that, for comparison with existing data, we take to be the mass of the charmed 
quark 1.275 GeV Fourier transform of the gluon propagator is 



^ ID 

A{p) = y ^—^^ (36) 



n=0 ' 



being 



So, we finally get 



yielding our final result 



2 rl)n+l -{n+^)n 



(^C?.G) = a, 0.2625 + 0.3304-^) (38) 



-C? • G) = (0.0039 ± 0.0003)GeT/'' (39) 

TT 

having taken as(m,-) — 0.34 ± 0.03. This is fully consistent with the bound 

(— G • G) < O.OOSGey* (40) 

TT 

and values currently known . But this is an evidence obtained directly from QCD of the existence of the gluon 

condensate that gives a firm ground to all the conclusions drawn assuming its existence. Let us point out that some 
recent analysis point toward an higher value for the gluon condensate [25j and this question is currently matter of 
debate. 

We have presented a proof of the observability of the pure Yang-Mills spectrum in current experiments. The 
theorem so proved permitted to obtain an effective theory to compute the width of the a resonance that in this way 
appears clearly as the ground state of Yang-Mills theory. The obtained width is in close agreement with experimental 
data for proper values of the strong couplig constant. Finally, we showed the existence of a gluon condensate that 
firmly estabilishes the identification of a already obtained in this way. 
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